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Dirac semimetals host bulk band-touching Dirac points and a surface Fermi loop. We develop 
a theory of superconducting Dirac semimetals. Establishing a relation between the Dirac points 
and the surface Fermi loop, we clarify how the nontrivial topology of Dirac semimetals affects their 
superconducting state. We note that the unique orbital texture of Dirac points and a structural phase 
transition of the crystal favor symmetry-protected topological superconductivity with a quartet of 
surface Majorana fermions. We suggest possible application of our theory to recently discovered 
superconducting states in Cd 3 As 2 . 


Dirac semimetals are three-dimensional (3D) materials 
that possess gapless (Dirac) points in the bulk Brillonin 
zone (BZ), whose low-energy excitations are effectively 
described as Dirac fermions. With time-reversal sym¬ 
metry (TRS) and inversion symmetry (IS) preserved, a 
pair of Dirac points is formed at the crossing of two dou¬ 
bly degenerate bands on a high-symmetry axis. The:v 
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are protected by discrete rotation {Cn) symmetry 
which prohibits band mixing to open a gap. Further¬ 
more, Dirac semimetals may host a surface Fermi loop 
(FL) [1-0. This contrasts sharply with a surface Fermi 
arc in Weyl semimetals Q because its topological origin is 
different. Several Dirac semimetals, including NasBi [1- 
ini and Cd 3 As 2 iiliilil , have been demonstrated 
experimentally and predicted theoretically [Slii. 

Superconducting phase transitions were reported re¬ 
cently in Cd 3 As 2 [23|-[25| and Au 2 Pb [ll], both of which 
support Dirac points protected by C 4 symmetry. Bulk 
Cd 3 As 2 exhibits superconductivity under high pressure 
(^8.5 GPa) [m accompanied by a structural phase tran¬ 
sition of the crystal [27| . In addition, point contact mea¬ 
surements of Cd 3 As 2 reportedly induce superconductiv¬ 
ity around the point contact region, where the tunneling 
conductance shows a zero-bias conductance peak [illii. 
Au 2 Pb also exhibits a superconducting phase transition 
after a structural phase transition [26| . 

In this letter we address the effect of the nontrivial 
topology, i.e., the Dirac points and FL, on the supercon¬ 
ducting properties. Topological materials are a promising 
platform to realize topological superconductors (TSCs) 
owing to the nontrivial topology of the wave function 
in normal states [isl-lisj. For instance, surface Dirac 
fermions may realize a TSC even for an s-wave pairing 
Also, the Fermi surface topology, which 
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state 

is the simplest topological structure in the normal state, 
directly affects the topological superconductivity of odd- 
parity superconductors [30|, [SJ. For the carrier-doped 
topological insulator, topological superconductivity has 
been anticipated for the surface [20 or the bulk [30. 

Here we present a general framework for studying su¬ 
perconductivity in Dirac semimetals. The key ingre¬ 
dients are symmetry-protected topological numbers in 
crystalline insulators and superconductors (iil-libj. In 


particular, we examine the C 4 topological invariant and 
the mirror Chern number, which ensure the existence of 
Dirac points and FLs, respectively, in Dirac semimet¬ 
als. First, we show that these two topological numbers 
are intrinsically related to each other, establishing a re¬ 
lation between Dirac points and surface FLs. Then, we 
elucidate how the nontrivial topology of Dirac semimet¬ 
als affects their superconducting state. We find that for 
a class of pairing symmetries, Dirac points and FLs in 
the normal state are inherited as bulk point nodes and 
surface Majorana fermions (MFs), respectively, in the su¬ 
perconducting state. 

By carefully examining the low-energy effective Hamil¬ 
tonian, we also reveal that doped Dirac semimetals favor 
an equal-spin odd-parity pairing rather than a conven¬ 
tional 5- wave one. The former pairing exhibits a distinct 
quartet of surface MFs stemming from the FL, though 
point nodes exist when the system retains C 4 rotation 
symmetry. If the C 4 symmetry is reduced to C 2 by a 
structural phase transition, the nodes disappear, and a 
full-gapped symmetry-protected TSC is realized. The 
FL-induced MFs are clearly distinguished from those in 
other TSCs [ 13,110 including superfluid ^He-A [ 10 , [50 
and Weyl superconductors [33|, I5ll454| . We finally suggest 
possible application of our theory to recently discovered 
superconducting states in Cd 3 As 2 and Au 2 Pb. 

Stability of Dirac points and surface Fermi loop — 
First, we provide a general argument on the topology of 
Dirac points. Our theory assumes TRS, IS, and uniaxial 
rotation symmetry, which are the most common sym¬ 
metries for Dirac semimetals. In the presence of TRS 
and IS, Kramer’s degeneracy exists at arbitrary k in the 
BZ, ensuring fourfold degeneracy when the conduction 
and valence bands are in contact. Such an accidental 
band crossing is generally not stable owing to band re¬ 
pulsion. However, if the band-touching point is on the 
high-symmetry axis, a Cn wmmetry can retain the band 
crossing as a Dirac point [ll-10 . Below we clarify the rel¬ 
evant topological structures. 

We focus on the C 4 symmetric Hamiltonian, 
C 4 H{k)C^^ = H{R^k), with R/^k = {ky, —k^, kz). Here, 
without losing generality, we have chosen the rotation 
axis as the kz axis. The commutation relations between 
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FIG. 1. (Color online) (a) Schematic picture of the Dirac 
points at (0, 0, =b/co) and the surface Fermi FL in the bulk BZ 
and the surface BZ, respectively, (b) On the /c^-axis, a bulk 
Dirac point with Q = — 1 in the left-hand side is decomposed 
into four chiral modes with different ap in the right-hand side. 
The bold line is doubly degenerate, (c) Schematic illustration 
of the surface Majorana quartet in the superconducting Dirac 
semimetal. 

the TRS operator T, IS operator P, and C4 are sum¬ 
marized as [T, P] = [T, C4] = [P, C4] = 0. As illus¬ 
trated in Fig. m^a), for a C4 symmetric tetragonal crys¬ 
tal, there are two C4 symmetry lines, FZ = ( 0 , 0 , kz) and 
MA = (tt, tt, /c^), with kz G [—tt, tt]. On these C4 symme¬ 
try lines, the Hamiltonian commutes with C4; thus, any 
energy band along the high-symmetry lines has a defi¬ 
nite eigenvalue of the C4 operator, = exp (p + |)] 
(p = 0 ,l, 2 , 3 ). 

The existence of Dirac points is ensured by the topo¬ 
logical invariant defined below. The system is generally 
gapped at the C4 symmetry points /c = F, Z, M, A, so one 
can count the number of bands below the Fermi level at 
these points. Denoting the number of such bands with 
the C4 eigenvalue at k as Np{k), we can introduce 
the topological number, Qp = Np{T) — Np{Z), for Dirac 
points on the FZ line, which we call the C4 topological in¬ 
variant. The Kramer’s degeneracy due to PT symmetry 
requires Qp = Qs-p (p = 0 , 1 ). Moreover, the sum rule 
= 0 holds for Dirac semimetals. Indeed, for a 
Dirac semimetal to become an insulator by a small C4 
breaking perturbation, the total number of bands below 
the Fermi level should be the same at all the symmetry 
points /c, which leads to the sum rule. From these two 
relations, we have 

Qo = —Qi = —Q2 = Qs = Q‘ (1) 

If Q 7^ 0, any band with ap has Q gapless points on the 
FZ lines, corresponding to the difference Np(r) — Np(Z), 
which eventually form Q Dirac points. In Fig. [IJb), we 
illustrate a Dirac point with Q = —1. Because band mix¬ 
ing between different C4 eigensectors is prohibited, the 
resultant Dirac points are stable as long as C4 symmetry 
is maintained. Similarly, we can also introduce the C4 


invariant for Dirac points on the MA line. 

Because of IS and the C2 subgroup for C4 symme¬ 
try, the system also has mirror reflection symmetry, 
MxyH(^kxiky^kz)M^y^ = Hikx’jky^ kz)-, with Mxy = 
C^P. In the mirror-invariant planes {kz = 0 or tt), the 
Hamiltonian is block-diagonal in the basis of the eigen¬ 
states of Mxy] thus, the mirror Chern number vx{kz) with 
kz = 0,7r is defined as h'x{kz) := l/ 27 r dkx dky p^{k)^ 
with A^{k) := T,E^<o^i^nik)\dk,u^{k)) where 

\u^{k)) is an eigenstate of H{k) in the mirror sector with 
the eigenvalue A = of Mxy and is the field strength 
of A^. Generalizing the relation between a band inver¬ 
sion and the Chern number in terms of eigenvalues of 
crystal symmetry (dH, , we obtain the following re¬ 

lation for the mirror Chern number: 

g^^A(O) ^ J^^P(P,A)(r) + Ar(p,A)(M)] ^2) 

p q 

where = exp {q = 0,1) is the eigenvalue of 

C'2, {k) is the number of occupied bands at k with a 

set of the C4 and Mxy eigenvalues A), and A/(g,A) (^) 
are those at Y with a set of C2 and Mxy eigenvalues 
(^g,A). [Y is the C2 symmetry point in Fig. [TJa)]. Note 
that occupied bands at the C4 (C2) symmetry points have 
a definite set of eigenvalues for C4 {C2) and Mxy because 
[(^4, Mxy] = [C2, Mxy] = 0 . We can also obtain a similar 
relation for z/A(7r) by replacing F, M, and Y with Z, A, 
and T, respectively, in Fig. [TJa). 

To see the close relationship between the C4 invari¬ 
ant and the mirror Chern number, consider a process 
in which a pair of stable Dirac points is created at F 
0 - Band inversion at F occurs in this process, so a 
Kramer’s pair of occupied bands, which have a set of 
eigenvalues (o^p/. A) and (^3-^/, —A) = A*), go above 

the Fermi level, and a Kramer’s pair of empty bands 
with the eigenvalues {apff, X) and (<a3-p",—A) go below 
it at F. As a result, N(^p A)(r) changes by AA'^p/ A)(r) = 
A7V(3-pq-A)(r) = -1, AiV^p^A)(^) = A7V(3-p",-A)(r) = 
1. To have a stable Dirac point, Qp should change at the 
same time, so p' ^ p" . Then, from Eq. ([ 2 ]), we find that 
this process induces a simultaneous change in the mirror 
Chern number Az/a( 0 )- 

^*^a( 0) = T] (p+^Ar(p^A)(r) mod 4 , 

p \ / 

= {p" — p') ^ 0 mod 4 . ( 3 ) 

Therefore, the creation of stable Dirac points is always 
accompanied by a net change in the mirror Chern num¬ 
ber. From the bulk-boundary correspondence, the resul¬ 
tant mirror Chern number ensures the existence of sur¬ 
face helical Dirac fermions, whose Fermi surfaces form 
FLs. 

Although the surface FL accompanies bulk Dirac 
points, it can be stable even when C4 symmetry is lost. 
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so the Dirac points have gaps. Indeed, unless the C2 sub¬ 
group is broken, the system maintains mirror reflection 
symmetry, which is sufficient to stabilize the surface FL. 
Therefore, the structural phase transition that breaks (^4 
to C2 retains the FL. 

Topology of superconducting Dirac semimetals — With 
a finite carrier density, Dirac semimetals have discon¬ 
nected bulk Fermi surfaces, each of which surrounds one 
of the band-touching Dirac points. See Fig. [T]^a). Now 
consider a superconducting state in Dirac semimetals. 
The system is described by the Bogoliubov-de Gennes 
(BdG) Hamiltonian, 


_ fH{k) - ^ A{k) \ 

HBdG[k) ^H*{-k)+n)’ 

where H{k) is the Hamiltonian for Dirac semimetals dis¬ 
cussed above, // is the chemical potential corresponding 
to the finite carrier density, and A{k) is the gap function. 
The BdG Hamiltonian supports particle-hole symmetry, 
CHBdG{k)C~^ = —iLBdG(~^): C = T^K^ with the Pauli 
matrix in the Nambu space and the conjugation op¬ 
erator K. Moreover, it may retain the symmetries of 
Dirac semimetals, depending on the symmetry property 
of the gap function. In particular, for a gap function 
with C4,/S.(k)C\ = c~^ A{R4^k) (r = 0 , • • • , 3 ), i?BdG(^) 
keeps C4 symmetry, C4.HBdG{k)C^^ = i^BdG(^4^), with 
C4 = diag[C4, e^Cl]. In addition, for a mirror-even or 
mirror-odd gap function that satisfies MxyA(k)Mly = 
riMA(kx^ky^—kz) with 77M = ± 1 , the system re¬ 

tains mirror reflection symmetry, MxyHBdG{k)M~y = 
HBdG{kx,ky, -kz), with M^y = diag[Ma,^, Gor- 

respondingly, we can introduce the C4-invariant and 
the mirror Ghern numbers Dx for iLBdG(^) [ 40 , 1 ^, in 
a manner similar to that used for those of H(k). The 
topological numbers Qp and z>a are responsible for the 
existence of bulk point nodes on the F Z line and surface 
MFs in the superconducting state, respectively. 

To evaluate these topological numbers, we employ the 
weak pairing assumption [SOj, [Sl], [SJ , i.e., that the super¬ 
conducting gap is much smaller than the Fermi energy. 
The gap function is reasonably negligible away from the 
Fermi surface, in which we can take A{k) 0 , leading 
to i^BdG(^) ^ didig[H{k) — la, —H*{—k) + //]. There¬ 
fore, at the symmetry points k = F, Z, M, H, we can 
relate the negative energy states of i^BdG(^) fo those 
of H(k). By taking into account the contribution from 
holes as well as electrons, the number Np of negative 
energy states with the eigenvalue ap is evaluated as 
Np{k) = Np(k)-^[N — Np^(k)], where N is the total num¬ 
ber of bands in H{k)^ ph = 3 — p + r mod 4 , and the 
first (second) term on the right-hand side comes from the 
electron (hole) contribution. From this equation, the (^4 


invariant in the superconducting state is obtained as 

g fo r = 0o.,=,». 

\2Qp otherwise. 

Similarly, the mirror Ghern number in the superconduct¬ 
ing state is calculated as the sum of the electron and hole 
mirror Ghern numbers, z>a = with \h = — 

As TRS in Dirac semimetals implies z/_a = —z^a, we have 


kx 


0 = 1 , 

2z^a = —1- 


(6) 


Relations dSD and have important physical con¬ 
sequences. (i) In the presence of C4 symmetry, any 
superconducting Dirac semimetal with a nontrivial r 
(r = 1 , 2 , 3 ) hosts point nodes as a remnant of Dirac 
points. Indeed, Eqs. o and ^ imply that at least a cou¬ 
ple of Qp are nonzero in this case. To open a point node 
gap, we need to break the C4 rotation symmetry, (ii) If 
the gap function is mirror-odd, the mirror Ghern num¬ 
ber of the superconducting Dirac semimetal is nonzero, 
resulting in double MFs. In Fig. [U^c), we illustrate how 
the double MFs are created. In general, the gap function 
mixes the surface FL of electrons with that of holes so as 
to open a gap for the FL. In the mirror-odd case, how¬ 
ever, mixing is prohibited on the mirror-invariant line in 
the surface BZ, so a pair of gapless points remains for 
each FL, forming double MFs. 

In addition to the double MFs on the mirror-invariant 
line, we also find that each FL in the mirror-odd super¬ 
conductor creates another pair of MFs on the kz axis 
in the surface BZ. By combining with C and Mxy, the 
BdG Hamiltonian for the surface FL on the kz axis, 
which we denote H^^Q{kz), has antiunit ary antisym- 
metry, CMHBdG{kz)C]^ = with Cm = 

iMxyC = iMxyTxK. Because = 1 in the mirror- 
odd superconductor, H^^Q{kz)iMxyTx is found to be real 
antisymmetric; thus, by using the Pfafhan, we can intro¬ 
duce the zero-dimensional topological invariant x{kz) = 
^gn{Pi[H^^Q{kz)iMxy'^x\}^ In the weak pairing case, 
x{kz) is evaluated as x{kz) = sgii{det[H^^{kz) — A^]}, 
where H^^{kz) — a is the Hamiltonian of the surface FL 
on the kz axis [^. Therefore, xi^z) has different signs 
inside and outside the FL, which implies that H^^Q{kz) 
should have zero-energy states near the points of intersec¬ 
tion between the FL and the kz axis. These zero-energy 
states form a pair of MFs on the kz axis. Gonsequently, 
we can conclude that each FL has a quartet of MFs, as 
shown in Fig. [U^c). The quartet of MFs can stay gap¬ 
less even when C4 symmetry is broken, as long as mirror 
symmetry is preserved. 

Low-energy analysis and application to CdsAs2. —For 
definiteness, we study the low-energy effective Hamil¬ 
tonian, which describes a class of Dirac semimetals in¬ 
cluding Cd3As2 and Au2Pb. Because bands in Dirac 
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semimetals are doubly degenerate owing to PT symme¬ 
try, band-touching Dirac points are minimally described 
by a 4 X 4 matrix Hamiltonian. Thus, in the minimal 
setup, we need orbital degrees of freedom in addition to 
spin degrees of freedom, which are given by the Pauli 
matrices cr^ and in the orbital (1,2) and spin (t^i) 
spaces, respectively. The form of the 4x4 Hamiltonian 
is uniquely determined by symmetry Hi- In particular, 
for P = Paz, the low-energy lattice Hamiltonian is given 
by 0 

H{k) ={M — txy{cos kx + cos ky) — tz cos kz}(TzSo 
P {rj sin kx)(TxSz - {rj sin ky)aySo 
P {/3 P^) sin kz (cos ky — cos kx)(TxSx 
- - ^){sinkzSinkxSinky)axSy, (7) 

with T = idoSyK and C 4 = . Here M, txy, 

tz, Tj, P, and 7 are material-dependent real constants. If 
tz > {M — 2txy) > 0, this model has a pair of Dirac 
points located nt k = ( 0 , 0 , Pko), with /cq > 0 defined by 
M = tz cosko -h 2txy We find z^±i( 0 ) = ±1 and \Q\ = 1 . 
Accordingly, the FL arises at the surface parallel to the 
kz axis, and the Dirac points are protected by C 4 . 

Near the Dirac points nt k = (0,0,±/co), the low- 
energy Hamiltonian takes the form of the Dirac Hamilto¬ 
nian, H{}Pj — Ptzk^ikz T ^0)^P z ^y^y^Op 

which exhibits nontrivial orbit-momentum locking. In 
Fig. [2fa), we show orbital textures in the kxky plane 
with kz = Pko in each spin sector, where the orientation 
of the orbit is tightly locked to the direction of the mo¬ 
mentum on the Fermi surface. Orbit-momentum locking 
critically affects the possible pairing symmetry in the su¬ 
perconducting state. Indeed, we can show that constant 
5 -wave pairing is inconsistent with the orbital texture. 
First, in such a static pairing state, electrons forming the 
Cooper pair have opposite momentum to each other, so 
they must belong to different Dirac points. Furthermore, 
as an s-wave pairing is spin-singlet, it must be formed 
between electrons in different spin sectors. However, for 
a Cooper pair between electrons in different Dirac points 
and different spin sectors, orbit-momentum locking re¬ 
quires a momentum-dependent orbital structure in the 
Cooper pairing, as illustrated in Fig. [2fa). Therefore, 
even if the pairing interaction favors an s-wave supercon¬ 
ducting state, the pairing function cannot be constant, 
suggesting suppression of the critical temperature. 

On the other hand, for a Cooper pair with parallel 
spins, orbit-momentum locking is consistent with a con¬ 
stant pairing function. Indeed, the orbital-singlet equal- 
spin pairing, A = Ao(c-f,iCt,2 - c^,2C4,.i) + iAo(ct,iCt,2 + 
C4.,2C4.,i)(= A||) is compatible with the orbital texture 
in Fig. Ha). Such an orbital-singlet Cooper pair is re¬ 
alized when the effective pairing interaction is domi¬ 
nated by an attractive interorbital interaction 'Hint = 
-2Vniri2, with = Es=t.4, 4 ,<T(®)cs,<T(a:) {V > 0 ) 
|3Q| . I62l - l64| . Although the actual pairing interaction is 



FIG. 2. (Color online) (a) Orbital textures on the Fermi sur¬ 
faces in the kxky-p\a,ne with kz = Pko in the = I (left) 
and = — I (right) sectors. Arrors indicate the direction of 
the orbit, {{(Tx), (c^))- A Cooper pair between electrons with 
opposite spin may realize both of a parallel orbit pair (1) and 
an anti-parallel one (2), depending on the momentum. On 
the other hand, a Cooper pair between electrons in the same 
spin state always have the anti-parallel orbit configuration (3). 
(b) Energy spectra at the (100) face, fi/tz = 0.5, M/tz — 4, 
txy/tz = 2, g/tz = I, P/tz = 2, 'g/tz = I, Ao/tz = O.I, 
A'oItz = 0.01, and mojtz =0.2. The distance between left 
{x — 0) and right {x — L) surfaces is L = 50. The red and 
green lines show the FL induced MFs. 


material-dependent, the above results imply that doped 
Dirac semimetals favor the latter gap function. Because 
AII is C 4 -symmetric with r = 2 and mirror-odd, i.e., 
C 4 A 11 C 4 = — A||, it realizes a symmetry-protected TSC 
with bulk point nodes and a surface MF quartet, as dis¬ 
cussed previously. In Fig. [2fb), we illustrate the quar¬ 
tet of MFs in this phase by numerically calculating the 
surface energy spectrum of the BdG Hamiltonian with 
Eq. ([7j) and A = Ay. Here we have also taken into 
account a symmetry-lowering effect from C 4 to C 2 by 
phenomenologically adding {mo sin kz)(TxSx to Eq. Q. 
As expected. Fig. Efb) proves the existence of MFs on 
the mirror-invariant line (FT) and the kz axis {TZ) and 
moreover shows a gap in the ZT direction due to the C 4 
breaking term [^. The obtained MFs in the mirror in¬ 
variant plane stay gapless even if the interaction effects 
are take into account 6 ^ . 

Finally, we discuss possible application of our theory 
to superconductivity in Au 2 Pb and Cd 3 As 2 . For Au 2 Pb, 
first-principle calculations show that the Fermi level of 
this material is inside the gap of the Dirac points ( 2 ^ : 
thus, no electron near the Dirac points contributes to the 
Cooper pairs. Hence, no TSC as discussed above is ex¬ 
pected in Au 2 Pb. On the other hand, the analysis above 
is applicable to the recently discovered superconductor 
Cd 3 As 2 . In Cd 3 As 2 , under high pressure, the structural 
phase transition occurs before the superconducting tran¬ 
sition. Together with the orbit-momentum locking dis¬ 
cussed above, the symmetry-lowering effect may stabilize 
the TSC phase by increasing the condensation energy, as 
the point nodes in the TSC phase are gapped when C 4 
is reduced to ( 72 - Therefore, it is likely that Cd 3 As 2 
realizes the TSC phase. The mirror-odd gap function 
of the TSC is detectable via anomalous Josephson ef- 
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fects 3^, @ with carefully fabricated junctions. 
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